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ABSTRACT 

We present the data analysis and the X-ray source counts for the first season of XMM- Newton ob- 
servations in the COSMOS field. The survey covers ~2 deg 2 within the region of sky bounded by 
9 ft 57.5 m < R.A. < I0 h 03. 5 m ; l d 27.5 m < DEC < 2 d 57.5 m with a total net integration time of 504 
ks. A maximum likelihood source detection was performed in the 0.5-2 keV, 2-4.5 keV and 4.5-10 
keV energy bands and 1390 point-like sources were detected in at least one band. Detailed Monte 
Carlo simulations were performed to fully test the source detection method and to derive the sky 
coverage to be used in the computation of the logN-logS relations. These relations have been then 
derived in the 0.5-2 keV, 2-10 keV and 5-10 keV energy bands, down to flux limits of 7.2xl0 -16 erg 
cm~ 2 s _1 , 4.0xl0 -15 erg cm~ 2 s _1 and 9.7xl0 -15 erg cm -2 s _1 , respectively. Thanks to the large 
number of sources detected in the COSMOS survey, the logN-logS curves are tightly constrained over 
a range of fluxes which were poorly covered by previous surveys, especially in the 2-10 and 5-10 keV 
bands. The 0.5-2 keV and 2-10 keV differential logN-logS were fitted with a broken power-law model 
which revealed a Euclidean slope (a ~2.5) at the bright end and a flatter slope (a ~1.5) at faint 
fluxes. In the 5-10 keV energy band a single power-law provides an acceptable fit to the observed 
source counts with a slope a ~2.4. A comparison with the results of previous surveys shows good 
agreement in all the energy bands under investigation in the overlapping flux range. We also notice a 
remarkable agreement between our logN-logS relations and the most recent model of the XRB. The 
slightly different normalizations observed in the source counts of COSMOS and previous surveys can 
be largely explained as a combination of low counting statistics and cosmic variance introduced by 
the large scale structure. 

Subject headings: cosmology: observations — cosmology: large scale strutcure of universe — cosmol- 
ogy: dark matter — galaxies: formation — galaxies: evolution — X-rays: surveys 



1. INTRODUCTION 

The source content of the X-ray sky has been in- 
vestigated over a broad range of fluxes and solid an- 
gles thanks to a large number of deep and wide surveys 
performed in the last few ye ars using ROSAT, Chan- 
dra and XMM-Newton (see iBrandt fc Hasingerl 120051 
for a review). Follow-up observations unambiguously 
indicate that Active Galactic Nuclei (AGN), many of 
which are obscured, dominate the global energy output 
recorded in the cosmic X-ray background. The impres- 
sive amount of X-ray and multi-wavelength data ob- 
tained to date have opened up the quantitative study 
of the demography and evolution of accretion driven 
Supermassive Black Holes (SMBHs: IM vaii et al.l I2000L 
lHasinger et~ai] l2005t lUeda et alj l2003t lLa Franca et all 
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2005). At present the two deepest X-ray surveys, the 
Chandra Deep Field North (CDFN; iBauer et al.ll2004h 
and Chandra Deep Field South (CDFS; iGiacconi et al.l 
120011 ). have extended the sensitivity by about two or- 
ders of magnit ude in all bands wi t h respect to pre- 
vious surveys (lHasinger et all 119931 : lUeda et all 119991: 
iGiommi et al.l 120001 ) . detecting a large number of faint 
X-ray sources. However, deep pencil beam surveys are 
limited by the area which can be covered to very faint 
fluxes (typically of the order of 0.1 deg 2 ) and suffer 
from significant field to field variance. In order to cope 
with such limitations, shallower surveys over larger ar- 
eas have been undertaken in the last few y ears with both 
Chan dra (e.g. the 9 deg 2 Bootes survey ( Murray et alJ 
2005), the Extended Groth strip EGS (jNandra et all 
20051) . the Extended Chandra Deep Field South E- 



CDFS, (iLehmer et all 120051: IVirani et alJlgOOll) and the 
Champ (jGreen et alJl2004t iKim et al.ll2004D )~ and XMM- 
Newt on (e.g. the HELLAS2X MM survey dFiore et all 
l200l . the XMM-Newton BSS dDella Ceca et al.ll20'ol 
and the ELAIS SI survey (jPuccetti et al.ll2006f) ). 

In this context the XMM-Newton wide field survey 
in the COSMOS field (IScoville et all 12001 . hereinafter 
XMM-COSMOS (|Hasinger et all 120071) . has been con- 
ceived and designed to maximize the sensitivity and sur- 
vey area product, and is expected to provide a major 
step forward toward a complete characterization of the 
physical properties of X-ray emitting SMBHs. A con- 
tiguous area of about 2 deg 2 will be covered by 25 in- 
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dividual pointings, repeated twice, for a total exposure 
time of about 60 ksec in each field. In the first ob- 
serving run obtained in A03 (phase A), the pointings 
were disposed on a 5x5 grid with the aimpoints shifted 
of 15' each other, so as to produce a contiguous pat- 
tern of coverage. In the second run, to be observed in 
A04 (phase B), the same pattern will be repeated with 
each pointing shifted by 1' with respect to phase A. The 
above described approach ensures a uniform and rela- 
tively deep coverage of more than 1 deg 2 in the cen- 
tral part of the field. When completed, XMM-COSMOS 
will provide an unprecedentedly large sample of about 
2000 X-ray sources with full multi-wavelength photo- 
metric coverage and a high level of spectroscopic com- 
pleteness. As a consequence, the XMM-COSMOS sur- 
vey is particularly well suited to address AGN evolution 
in the context of the Large Scale Structure in which they 
reside. More specifically, it will be possible to investi- 
gate if obscured AGN are biased tracers of the cosmic 
web and whether their space density rises in the proxim- 



ity o f galaxy clusters (IHenrv fc Briell 1991t ICappi et al 
20011 : iGilli et all [20031: 1 Johnson et alj l2003t lYang et al 



2003t ICappelluti et al.ll2005t iRuderman fc Ebeling||2005l: 
Mivaii et all 120071 : lYang et all 120061 : ICappelluti et all 
2006). 

The X-ray reduction of phase A data along with a de- 
tailed analysis of the source counts in different energy 
bands are presented in this paper which is organized as 
follows. In Section 2 the data reduction procedure and 
the relative astrometric corrections are described. In Sec- 
tion 3 the source detection algorithms and technique are 
discussed. Monte Carlo simulations are presented in Sec- 
tion 4. The logN-logS relations and the analysis of the 
contribution of the XMM-COSMOS sources to the X- 
ray background are discussed in Section 5. The study 
of sample variance is presented in Section 6 and a sum- 
mary of the work is reported in Section 7. The strategy 
and the log of the observati o ns of XMM-COSMOS are 
presented by IHasinger et alj (120071). the optic al identifi- 
cations of X-rav sources bv lBrusa et al. (2007), th e ana l- 
ysis of groups and clusters by iFinoguenov et all (|2007h , 
th e spectral analysis of a subsample of bright sources 
by iMainieri et al.l (|2007 ) and the clu s tering of X-ray ex- 
tragalactic sources by iMiyaii et al.l (|2007h . Through- 
out the paper the conco rdance WMAP ACDM cosmol- 
ogy (jSpergel et al.ll2003h is adopted with H o =70 km s _1 
Mpc" 1 , fi A =0.7 and ft m =0.3 

2. EPIC DATA CLEANING 

The EPIC data were processed using the XMM- 
Newton Standard Analysis System (hereinafter SAS) 
version 6.5.0. The Observational Data Files (ODF, "raw 
data") of each of the 25 observations, were calibrated us- 
ing the SAS tools epchain and emchain with the most 
recent calibration data files. Events in bad columns, bad 
pixels and close to the chip gaps were excluded. 
Both the EPIC PN and MOS event files were searched for 
high particle background intervals. The distribution of 
the background counts binned in 100 s intervals was ob- 
tained in the 12-14 keV band for the PN and in the 10-12 
keV band for the MOS, which are dominated by particle 
background, and then fitted with a gaussian model. All 
time intervals with background count rate higher than 3cr 
above the average best fit value were discarded. In Fig. 




counts/100 s 

Fig. 1. — The background counts distribution in the PN ob- 
servation of Field 6. The solid line represents the best gaussian 
fit to the distribution. The continuous vertical line represents the 
adopted 3<r cut above which the corresponding time intervals have 
been discarded. 



Q] an example of the application of this method to Field 
#6 is shown. Once the high energy flares were removed, 
the 0.3-10 keV background counts distribution was pro- 
cessed, with the same 3tr clipping method, in order to re- 
move times during which low energy particle flares were 
important. These flares are not easily detected in the 
12-14 keV band. As a result of this selection process the 
average time lost due to particle flares was <20% and 
2 obs ervations were completely lost (sec Ftasinger et al. 
12001 . 

An important feature observable in the background spec- 
trum of both MOS and PN CCDs is the Al Ka (1.48 
keV) fluorescent emission. In the PN background two 
strong Cu lines are also present at ~ 7.4 keV and ~ 
8.0 keV. Since these emission lines could affect the sci- 
entific results, the 7.2-7.6 keV and 7.8-8.2 keV energy 
bands in the PN and the 1.45-1.54 keV band (in PN 
and MOS) were excluded from the detectors events. Im- 
ages were then created in the 0.5-2 keV, 2-4.5 keV and 
4.5-10 keV energy bands with a pixel size of 4 arcsec. 
Single and double events were used to construct the PN 
images, while MOS images were created using all valid 
event patterns. Out-of-Time (OOT) events appear when 
a photon hits the CCD during the read-out process in 
the IMAGING mode. The result is that the x position 
of the event/photon is known, while the y position is 
unknown due to the readout and shifting of the charges 
at this time. For this reason artificial OOT event files 
were created. A new y coordinate is simulated by ran- 
domly shifting the event along the readout axis and per- 
forming the gain and CTI (charge transfer inefficiency) 
correction afterwards. For the PN, in full frame mode 
the OOT events constitute about 6.3% of the observ- 
ing time. Those files were filtered in the same way as 
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Fig. 2.— Shift between the PN+MOS mosaic and the MEGA- 
CAM catalog for each pointing of the XMM-Newton COSMOS 
field. 



the event files and the produced images were subtracted 
from the event images. Images were then added in order 
to obtain PN+MOS mosaics. For each instrument and 
for each observation, spectrally weighted exposure maps 
were created using the SAS task eexpmap, assuming a 
power law model with photon index T=2.0 in the 0.5-2 
keV band and T=1.7 in the 2-4.5 and 4.5-10 keV bands. 

2.1. Astrometry correction 

In order to correct the astrometry of our XMM- 
Newton observations for each pointing and for each in- 
strument, the produced X-ray source list (see next sec- 
tion) was comp ared with the MEGACA M catalog of the 
COSMOS field (|Mc Cracken et al.ll200l including all the 
sources with I magnitudes in the range 18-23. In order to 
find the shift between the two catalogues, an optical- X- 
ray positional correlation was computed using the likeli- 
hood algorithm included in the SAS task eposcorr. This 
task uses in a purely statistical way all possible coun- 
terparts of an X-ray source in the field to determine 
the most likely coordinate displacement. This method is 
independent of actual spectroscopic identifications, but 
all post-facto checks using, for example, secure spectro- 
scopic identifications, have demonstrated its reliability 
and accuracy. Using the magnitude range mentioned 
above, systematic effects introduced by bright stars and 
faint background objects are minimized. In the majority 
of the observations the shift between the three cameras 
turned out to be < 1" (i.e. much smaller than the pixel 
size of the images used here ). Since the shift between 
the EPIC cameras is negligible, a correlation between the 
joint MOS+PN source list and the optical catalog was 
calculated to derive the astrometric correction. For the 
23 pointings presented here, the shifts between the opti- 
cal and X-ray catalog are never larger than 3" , with an 
average shift of the order of Aa ~1.4" and AS —0.17". 
The average displacement in the two coordinates between 
the PN+MOS mosaic X-ray positions and the MEGA- 
CAM catalog sources for each pointing of the XMM- 
COSMOS field is shown in FigH The appropriate offset 
was applied to the event file of each pointing and im- 
ages and exposure maps were then reproduced with the 
corrected astrometry. 



3. EPIC SOURCE DETECTION 

3.1. Background modeling 

In order to perform the source detection a sophisti- 
cated background modeling has been developed. In X- 
ray observations the background is mainly due to two 
components, one generated by undetected faint sources 
contributing to the cosmic X-ray background and one 
arising from soft protons trapped by the terrestrial mag- 
netic field. For this reason two background templates 
were computed for each instrument and for each point- 
ing, o ne for the sky (vignetted) background (|Lumb et al.l 
2002) and one for instrumental and particles background 
(unvignetted) . To calculate the normalizations of each 
template of every po inting, we first perform ed a wavelet 
source detection (see lFinoguenov et al.ll2007l ) without so- 
phisticated background subtraction, then we excised the 
areas of the detector where a significant signal due to 
sources was detected. The residual area is split into two 
parts depending on the value of the effective exposure 
(i.e. higher and lower than the median value). Using 
the two templates, we calculate the coefficients of a sys- 
tem of two linear equations from which we obtain the 
normalizations of both: 
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where A and B are the normalization factors, M„' 2 the 
vignetted templates in the region with effective exposure 
higher and lower than the median , M^ v the unvignetted 
templates and C 1,2 are the background counts in the two 
regions. The region with effective exposure lower than 
the median (i.e. high vignetting, >7' off-axis) is domi- 
nated by the instrumental background, while the region 
with higher effective exposure is dominated by the sky 
background. Therefore with this method we have the 
advantage of better fitting the two components of the 
background. The standard method for estimating the 
background, based on the spline functions used in the 
XMM-Afewfon pipelines, returned in our case significant 
residuals. The excellent result of this technique can be 
seen in the signal-to-noise (SNR) map in Fig [3j despite 
the significant variations in exposure time and average 
background level from pointing to pointing, a rather ho- 
mogeneous signal-to-noise ratio is achieved across the 
whole mosaic. It is worth noting that also pixels with 
negative values are shown in the map; these are located 
where the background model is higher than the measured 
background. 

3.2. Maximum likelihood detection 

In each pointing the source detection was conducted 
on the combined images of the different instruments in 
the three energy bands mentioned above using the SAS 
tasks eboxdetect and emldetect. As a first step, the sliding 
cell detection algorithm eboxdetect was run on the im- 
ages in the three energy bands. In this procedure source 
counts were collected in cells of 5 x 5 pixels adopting a low 
threshold in the detection likelihood (i.e. likemin=4). 
The source list produced by eboxdetect was then used as 
input for emldetect. For all the sources detected with the 
sliding cell method this task performs a maximum likeli- 
hood PSF fit. In this way refined positions and fluxes for 




N. Cappelluti et al. 




Fig. 3. — Signal— to-noise ratio map in the 0.5—2 keV band of the XMM-Newton raster scan in the COSMOS field. The stretch of the 
color map corresponds to [-0.1<S/N<1 per pixel]. The scale has been chosen to enhance the SNR contrasts. If S is the raw (sources + 
background) image and B is the model background image, then the SNR map is obtained by SNR = . The image was smoothed with 
a gaussian filter with <r=2 pixels. Negative values are places where the measured background is smaller than the model background. 



the sources were determined. Due to the particular pat- 
tern of our observations fsee lHasinger et aT]|2007f ). the 
same source could be detected in up to 4 different point- 
ings. For this reason both eboxdetect and emldetect were 
run in raster mode. The source parameters (position and 
flux) were fitted simultaneously on all the observations 
where the source is observable, taking into account the 
PSF at the source position in each pointing. As likeli- 
hood threshold for the detection, we adopted the value 
det_ml=6. This parameter is related to the probability 
of a random Poissonian fluctuation having caused the 
observed source counts: 

det-ml = —InPrandom (3) 

In principle, the expected number of spurious sources 
could be estimated as the product of the probability for 
a random Poisson fluctuation exceeding the likelihood 
threshold times the number of statistically independent 



trials, N trial- For a simple box detection algorithm N tr iaJ 
would be approximately given by the number of indepen- 
dent source detection cells across the field of view. For 
the complex multi-stage source detection algorithm, like 
the one applied here, N tr iai cannot be calculated ana- 
lytically, but has to be estimated through Monte Carlo 
simulations. These simulations, which are discussed in 
Section 4, return a number of spurious sources of ^2% 
at the likelihood level chosen. All the sources were fit- 
ted with a PSF template convolv ed with a beta model 
l|Cavaliere &: Fusco-Femianolll976l) . Sources which have 
a core radius significantly larger than the PSF are flagged 
as extended {ext parameter > 0). 

A total of 1307, 735 and 187 X-ray sources were de- 
tected in the three bands. Of these sources, twenty-six 
were classified as extended. The analysis of the X-ray ex- 
tended sources in the COSMOS field is beyond the scope 
of this work; these sources are extensively discussed by 
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TABLE 1 

Summary of source detection 



Energy Band 


r a 


ECF b 




N sources' 1 


Single detections 




keV 




cts s — 1 /10 — 11 erg cm -2 s _1 


erg cm -2 s — 1 


All 


Point-like 




0.5-2 


2.0 


10.45 


7.2xl0~ 16 


1307 


1281 


661 


2-4.5 


1.7 


1.52 


4.0X10" 15 


735 


724 


89 


4.5-10 


1.7 


1.21 


9.7xl0" 15 


187 


186 


3 



a r is the average spectral index assumed in each band. 

b The energy conversion factor. Note that the ECF values in the second and third rows are the conversion factors 
from flux in the 2—10 keV and 5—10 keV bands to count rate in the 2-4.5 keV and 4.5-10 keV bands, respectively. 
The ECFs are computed by assuming as a mean spectrum an absorbed power-law with N£f=2.6xl0 20 cm -2 
and spectral index T=2.0 in the 0.5-2 keV band and T=1.7 in the 2-4.5 keV and 4.5-10 keV bands. 
c The flux of the faintest source. 

d The total number of sources detected for the entire sample and the point-like sample only. 
e the number of sources detected only in one band. 



lFinoguenovet~aT1 (|200l . A total of 1281, 724 and 186 
point-like sources were detected in the three bands down 
to limiting fluxes of 7.2xl0~ 16 erg cm -2 s _1 , 4.7xl0~ 15 
erg cm -2 s _1 and 9.7xl0 -15 erg cm -2 s _1 respectively. 
The minimum number of net counts for the detected 
sources is ~21, 17 and 27 in the three bands, respectively. 
A total of 1390 independent point-like sources have been 
detected by summing the number of sources detected in 
each band but not in any softer energy band. The num- 
ber of sources detected only in the 0.5-2 keV, 2-4.5 keV, 
4.5-10 keV bands are 661, 89 and 3, respectively. 

From the count rates in the 0.5-2 keV, 2-4.5 keV and 
4.5-10 keV bands the fluxes were obtained in the 0.5- 
2 keV, 2-10 keV and 5-10 keV bands respectively using 
the energy conversion factors (ECF) listed in Table[IJ to- 
gether with a summary of the source detection. The ECF 
values have been computed using the most recent EPIC 
response matrices in the corresponding energy ranges. 
As a model, we assumed power-law spectra with N# 
=2.6 xlO 20 cm" 2 , (corresponding to the average value of 
N H o ver the whole COSMOS field (jDickev fe Lockmanl 
1990)) and the same spectral indices used to compute the 
exposure maps without considering any intrinsic absorp- 
tion (see Table [l| . It is worth noting that the spectral 
indices and the absorptions of the individual sources can 
be significantly diff erent from the average values assumed 
here. In particular Mai nieri et alj (|2007f ) found that the 
spectral indices T of the XMM- COSMOS sources are in 
the range 1.5+2.5, in the CDFS lTozzi et alJ (|2006f ) mea- 
sured an average photon ind ex < T >~1.75 a nd sim- 
ilar values were obtained by iKim et ail (|2004h in the 
CHAMP survey. The mean spectrum assumed here is 
therefore consistent with the values measured up to now. 
By changing the spectral index of ±0.3 the ECFs change 
of 2%, 12% and 4% in the 0.5-2 keV, 2-4.5 keV and 
4.5-10 keV, respectively. 

4. MONTE CARLO SIMULATIONS 

In order to properly estimate the source detection effi- 
ciency and biases, de tailed Monte Carlo simulations were 
perfo rmed (see, e.g.. lHasinger et aT1ll993l ; [Loaring et al.1 
2005). Twenty series of 23 XMM-Newton images were 
created with the same pattern, exposure maps and back- 
ground levels as the real data. The PSF of the simu- 
lated sources was constructed from the templates avail- 
able in the XMM-Newton calibration database. The 
sources were randomly placed in the field of view ac- 
cording to a standard 0.5-2 keV logN-logS distribution 




Fig. 4. — The cumulative probability to detect a true 
source with det_ml>6 in a circle of a given radius in the 0.5-2 
(continuous line), 2—4.5 (dashed line) and 4.5-10 keV (dashed — 
dotted line) energy bands. The 68%, 90% and 99% levels are plot- 
ted as horizontal lines. 



(|Hasinger et all 12001 . This was then converted to a 
2-4.5 keV and 4.5-10 keV logN-logS assuming that all 
the sources have the same intrinsic spectrum (a power- 
law with spectral index T = 1.7). We then applied, to 
the simulated fields, the same source detection proce- 
dure used in the real data. ISchmitt fe Maccacarol (|1986f ) 
showed that with the threshold adopted here for source 
detection, which corresponds roughly to the Gaussian 
4.5-5(7, the distortion of the slope of the logN-logS due 
to Poissonian noise is <3% for a wide range of slopes. 
Therefore, the uncertainties introduced by using a single 
logN-logS as base for the simulations are negligible. A 
total of 30626, 13579 and 3172 simulated sources were de- 
tected in the 0.5-2 keV, 2-4.5 keV and 4.5-10 keV bands 
down to the same limiting fluxes of the observations. For 
every possible pair of input-output sources we computed 
the quantity 



R 2 = 1 ''' ~ ' r,i 



y-yo 



0~ ii 



S-S 
OS 



(4) 



where x, y and S are the position and flux of the detected 
source and Xq, and Sq are the corresponding values for 
all the simulated sources We then flag as the most likely 
associations those with the minimum value of R 2 . The 
distribution of the positional offsets is plotted in Fig. [4] 
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for each energy band analyzed. 
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Fig. 5. — The ratio S otl t/Si„ as a function of the output de- 
tected flux in the 0.5-2 keV, 2-4.5 keV and 4.5-10 keV bands, 
respectivelyTop middle and bottom panel. 



from their input counterpart are classified as spurious. 
These account for 2.7%, 0.5%, 0.6% of the total num- 
ber of sources in the 0.5-2 keV, 2-4.5 keV and 4.5-10 
keV bands, respectively. Source confusion occurs when 
two or more sources fall in a single resolution element of 
the detector and result as a single detected source with 
an amplified flux. In order to determine the influence of 
the source confusion w e adopted the method described in 
lHasinger et alj (|1998l ). We define as "confused" sources 
those for which S out /(S in + 3 * a out ) > 1.5 (where a out 
is the la error on the output flux). The fraction of 
"confused" sources is 0.8%, 0.15% and <0.1% in the 
0.5-2 keV, 2-4.5 keV and 4.5-10 keV bands, respectively. 
The photometry was also tested; the ratio of output to 
input fluxes in the simulation is plotted in Fig|5] 

At bright fluxes this ratio is consistent with one, while 
at fainter fluxes the distribution of S ou t/Si n becomes 
wider, mainly because of increasing errors, and skewed 
toward values greater than one. This skewness of the 
distribution can be explained mainly by t wo effects: 
a) so urce confusion and b) Eddington Bias (Eddington 
fl940h . While source confusion, as defined above, affects 
only a small fraction of the sources, the Eddington bias 
results in a systematic upward offset of the detected flux. 
The magnitude of this effect depends on the shape of the 
logN-logS distribution and the statistical error on the 
measured flux. Since there are many more faint than 
bright sources, uncertainties in the measured flux will 
result in more sources being up-scattered than down- 
scattered. Together with this, the fact that in the 4.5-10 
keV band we are sampling a flux region in which the 
logN-logS is steeper than in the other bands (see Section 
5), explains why such an effect is more evident in the 
4.5-10 keV band. 

Besides assessing the reliability of our source detec- 
tion procedure, one of the aims of these simulations is 
to provide a precise estimation of the completeness func- 
tion of our survey, known also as sky coverage. We con- 
structed our sky-coverage (f2) vs. flux relation by di- 
viding the number of detected sources by the number of 
input sources as a function of the flux and rescaling it 
to the sky simulated area. Having analyzed the simula- 
tions with the same procedure adopted for the analysis of 
the data, this method ensures that when computing the 
source counts distribution (see next section) all the ob- 
servational biases are taken into account and corrected. 
The vs. flux relation relative to the 0.5-2 keV, 2-10 
keV and 5-10 keV bands is plotted in fig. [5] The total 
sky area is 2.03 deg 2 and it is completely observable down 
to fluxes of ~0.3, 1.3 and 2xl0 -14 erg cm~ 2 s _1 in the 
three bands, respectively. The sky coverage drops to at 
limiting fluxes of ^7xl0 -16 erg cm' 2 s _1 , ~4xl0~ 15 



1 &nd 



^9x10 15 erg cm 2 s~ 



erg 

\in the 0.5-2 keV, 



2-10 keV and 5-10 keV bands, respectively. 



We find that 68% of the sources are detected within 
2.1", 1.3" and 0.8" in the 0.5-2 keV, 2-4.5 keV and 
4.5-10 keV bands, respectively. Since the detection soft- 
ware fits the position of the source using the information 
available for the three bands together, we expect to be 
able to detect sources with an accuracy of the order of, 
or somewhat be t ter th an, that shown in Fig. 2J As 
m lLoaring et ail (|2005h . we then define a cut-off radius 
r cut of 6" . Sources with a displacement larger than r cut 



5. SOURCE COUNTS 

Once the sky coverage is known, the cumulative source 
number counts can be easily computed using the follow- 
ing equation: 

Ns 1 

A — 1 1 
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Fig. 6. — The sky coverage vs. flux relation in the 0.5-2, 2-10 
and 5-10 keV energy bands (respectively continuous, dashed and 
dash-dotted line), resulting from the simulations described in the 
text. 
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TABLE 2 

Cumulative number counts 



log(S) a 


n h 


N(>S) C 


erg cm 2 s 1 


deg 2 


deg - 2 


0.5-2 keV 


-13.0 


2.03 


4.5±1.5 


-13.5 


2.03 


18.8±3.1 


-14.0 


2.03 


105.2±7.0 


-14.5 


2.03 


327.0±12.7 


-15.0 


0.58 


790±23.3 


-15.1 


0.12 


931±53.0 


2-10 keV 


-13.0 


2.03 


8.6±2.0 


-13.5 


2.03 


57.0±5.3 


-14.0 


1.40 


258.9±11.6 


-14.3 


0.13 


600.1±34.2 


5-10 keV 


-13.5 


2.03 


21.3±3.2 


-14.0 


0.35 


lllill.O 


a Flux 






b Sky coverage 







b Cumulative source counts 

where N$ is the total number of detected sources in the 
field with fluxes greater than S and f2j is the sky coverage 
associated with the flux of the \ th source. The variance 
of the source number counts is therefore defined as: 



N s 

E 



(6) 



Source number counts are reported in Table O The 
cumulative number counts, normalized to the Euclidean 
slope (multiplied by S 1,5 ), are shown in Figures [3 [8] and 

El 

for the 0.5-2 keV, 2-10 keV and 5-10 keV energy 
ranges, respectively. With such a representation, the de- 
viations from the Euclidean slope are clearly evident as 
well as the flattening of the counts towards faint fluxes. 
Source counts are compared with the findings of other 
deep and shallow surveys collected from the literature. 
The plotted reference results were selected in order to 
sample a flux range as wide as possible and at the same 
time to keep the plots as clear as possible. As discussed 



Fig. 7.— The 0.5-2 keV logN-logS of the XMM-COSMOS 
(red dots) s ources compared wit h the ROSAT medium sensitiv- 
ity survey (Hasinger ct al. 1993) (blue do t dashed line), com - 
bined ROSAT, XMM-Newton , Chandra lIHasinger et alj |2005|) 
(gre en dashed line), the Chandra deep field south la error 
tie jRpsati ct al. 2002) (magent a continuous line ), the Chan- 
dra deep field north la error tie (Bauer ct al. 2004) (pink dot — 
dashed line), the 100 k s of the XMM— Newton Lockman 
hole lIHasinger et all [20011 (cyan circles), the HELLAS2XMM 
(Baldi ct al. 2002]) (black pentagons) and the extended CDFS 
IjLehmer et al.l 120051 ) (black continuous line) surveys. The over- 
layed bla ck-dashed line represents the logN- logS predicted by the 
model of IGilli. Comastri fc Hasinger I (2006h . The source number 
counts are plotted scaled by S 1 ' 5 in order to highlight the deviation 
from the Euclidean behavior. 



in the previous section, the sky coverage f2 was derived 
from realistic Monte Carlo simulations and therefore no 
further correction for the Eddington Bias is required. 

In order to parameterize our relations, we performed 
a maximum likelihood fit to the unbinned differential 
counts. We assumed a broken power-law model for the 
0.5-2 keV and 2-10 keV bands: 



f<t\ dN 



AS- ai S>S b 
BS- a2 S<S b , 



(7) 



where A — B S^ 1- " 2 is the normalization, ot\ is 
the bright end slope, a-i the faint end slope, Sb the 
break flux, and S the flux in units of 10~ 14 erg cm~ 2 
s _1 . Notice that using the maximum likelihood method, 
the fit is not dependent on the data binning and 
therefore we can make full use of the whole dataset. 
Moreover, the normalization A is not a parameter of 
the fit, but it is obtained by imposing the condition 
that the number of expected sources from the best 
fit model is equal to the total observed number. In 
the 0.5-2 keV energy band the best fit parameters are 

+0.28 xl0 -14 



ai=2.60±^|, a 2 =1.65±0.05, 56=1.55 



0.24 er S 
and A=123. Translating this value of the 
normalization to that for the cumulative distribution 
at 2xl0~ 15 erg cm~ 2 s _1 , which is usually used in 
the literature for Chandra surveys, we obtain Ai$ ^450 
which is fully consistent with most of previous works 
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Fig. 8.— The 2-10 keV logN-logS of the XMM-COSMOS 
(red dots) sources compared with Chandra, XMM-New ton and 
ASCA (blue dash ed line) IMoretti et al.l l2003h , HEL- 
LAS BeppoSAX IIGiommi et al] |200*oT) (black hexagons) 
the Chandra deep field south lcr error tie l|Rosati et al.l 
120021 ) (magenta continuous lin e), the HELLAS2XMM 
(green pentago ns) (Baldi ct al. 2002), the E LAIS SI (blue stars ) 
i'puccetti eFaTI 12009) . the extended CDFS (ILehmer et al.l 120051) 
lcr error tie (black continu ous line) and the 100 ks of the Lockman 
hole (cian open circles) JHasingcr ct al. 2001). The overlayed 
bl ack-dashed line represents the logN -logS predicted by the model 
of Gilli, Comastri & Hasingcr (2006). The source number counts 
are plotted scaled by S 1 ' 5 in order to highlight the deviation from 
the Euclidean behavior. 
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Fig. 9.— The 5-10 keV logN-logS of the XM M-COSMOS 
(red dots) sources compared with the HELLAS2XMM f Baldi et all 
12003) (green pentago ns), the Chandra deep field south 1(7 error 
tie feo sati ct al. 2002) (magenta continu ous line) , the HELLAS- 
BeppoSAX data from Fiorc et al. (2 001|) (bl ack hexagons), the 
ELAIS SI (blue stars ) ifPuccetti et alj 120061) and the 100 ks of 
the Lockman Hole (cyan open circles) (Hasingcr ct al. 2001). The 
overlayed bl ack-dashed line represents th e logN -logS predicted by 
the model of lGilli. Comastri fc Hasinger I ll200a ). The source num- 
ber counts are plotted scaled by S 1,5 in order to highlight deviation 
from the Euclidean behavior. 



where a fit result is presented (|Hasinger et alj 119931: 
Mushotzkv et al.H2qO0t lHasinger et al.|l200lUBaldi et al ' 
2002; iRosati et all 120021 : iBauer et~aU l2004t iKim et al 



2004t lHasinger et al.l 120051 : iKenter et al.ll2005h . but sig- 
nifica ntly lower t han t hat found in the CLASXS sur- 
vey (|Yang et alj 120041 ). In the 2-10 keV band the 
best fit parameters are ai=2.43±0.10, a2=1.59±0.33, 
SW-02±°;^ xl0~ 14 erg cm- 2 s^ 1 and ,4=266. The 
latest value translates into A 15 ~1250. Also in this band, 
our results are in agreement with previous surveys within 
lcr, with the exception of the CLASXS su rvey which is 
~30% higher in this band (|Yang et al.ll2004l ). In the 5-10 
keV energy bands, where the differential counts do not 
show any evidence for a break in the sampled flux range, 
we assumed a single power-law model of the form: 



H{S) = ~dJ 



(8) 



for which the best fit parameters are found to be 
,4=102 and ai =2.36±0.1. 

In the soft 0.5-2 keV (Fig. 8) energy range a visual 
inspection of the various datasets suggests a remarkably 
good agreement between XMM-COSMOS and literature 
data 7 . 

7 In particular, it is worthwhile t o notice the good agre ement 
between the XMM-COSMOS and the lHasinger et alj (2005) logN- 
logS, which has been used as input in our simulations. This good 
agreement can be considered as an a posteriori support of the re- 
liability of those results from the simulations which depend on the 



In the 2-10 keV band the XMM-COSMOS counts 
bridge the gap between deep field observations 
(|Rosati et al.l 120021) and sh allower large area Bep- 
poSA X (iGiommi et all I2000D and XMM-Newton sur- 
veys (iBaldi et al.l l2002f ). At relatively bright fluxes 
(> 10" 14 erg cm" 2 s" 1 ) the XMM-COSMOS logN-logS 
nicely matches previous measurements, though provid- 
ing a much more robust estimate of the source counts 
thanks to the much smaller statistical errors. 
A major step forward in the determination of X-ray 
source counts is achieved in the 5-10 keV band, where 
the previously existing data from different surveys show 
very significant differences. Thanks to XMM-COSMOS, 
a solid measure of the hard-X logN-logS in the flux in- 
terval 10~ 14 - 10 -13 erg cm -2 s _1 is obtained for the 
first time. From Fig. [3] we notice that the normaliza- 
tion of the XMM-COSMOS logN-logS is slightly higher 
than (~10%), although consistent at la with, that mea- 
sured by Chandra while, ELAIS SI (|Puccetti et al.ll2006|) 
source counts are 30% lower. However, in the overlapping 
flux range the latter is characterized by large errors due 
to the small number of relatively bright sources in the 
Chandra deep fields. Interestingly enough, the XMM- 
COSMOS counts match nicely, with smaller errors, those 
of the wide area HELLAS2XMM survey (Baldi et al. 
2002), while the p ioneering measurements of BeppoSAX 
(|Fiore et alj 120011 ) are systematically higher than the 
counts from XMM-COSMOS. 

assumed input logN-logS. 
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5.1. Resolved fraction of the X-ray background 

One of the main aims of XMM-COSMOS with its large 
and medium-deep coverage, is to provide a solid census 
of the X-ray source population to be compared with ob- 
servations and models of AGN ev olution. According to 
recent synthesis models (see e.g. ; IComastri et al.lH995t 
iGilli. Comastri fc Hasinger ||2006t IWorslev et al.ll2005fl a 
high fraction of heavily obscured AGN is necessary to 
explain the spectral shape and the intensity of the X-ray 
background (XRB). We examine therefore which fraction 
of XRB is resolved into discrete sources in our survey. 

As a first test, we computed the flux which XMM- 
COSMOS itself resolves into discrete sources by sum- 
ming their fluxes weighted on the sky coverage in the 
0. 5-2 keV, 2-10 keV a nd 5-10 keV energy bands. As 
m IWorslev et all (|2005t) we used as reference value of 
the normalization at 1 keV of the XRB spectrum that of 
iDe Luca fc Molendil (|200l which assumes that the spec- 
tral shape in the 1-10 keV band is a power-law with spec- 
tral index T=1.4 and a normalization at IkeV of 11.6 keV 
cm" 2 s" 1 keV" 1 . The latter value corresponds to a flux 
of 0.80, 2.31 and 1.27X10" 11 erg cm" 2 s" 1 deg" 2 in the 
0.5-2 keV, 2-10 keV and 5-10 keV energy bands, respec- 
tively. In the 0.5-2 keV band we measure a contribution 
of the sources to the XRB which corresponds to a normal- 
ization at 1 keV of 0.49±0.08xl0 _11 erg cm" 2 s" 1 deg -2 . 
The corresponding values in the 2-10 keV and 5-10 keV 
bands are 0.92±0.22 x 10 -11 erg cm" 2 s" 1 deg -2 and 
0.28±0.15xl0~ n erg cm" 2 s" 1 deg" 2 . Therefore XMM- 
COSMOS resolves by itself -65%, -40% and -22% of 
the XRB into discrete sources in the 0.5-2 keV, 2-10 keV 
and 5-10 keV energy bands, respe ctively. It is worth 
notici ng that the flux measured by IDe Luca fc Molendil 
(|2004f ) is the highest measu red in literature in the 1-10 
keV energy range (see e.g. IGilli. Comastri fc Hasingerl 
I2006L for a complete collection). It is also worth notic- 
ing that we computed the fraction of resolved XRB by 
assuming that all the sources have the same spectrum. 
Therefore, in our estimate, effects due to the broad ab- 
sorption and spectral index distributions of AGNs are 
not included. 

In Figs. El [5] and M we compared our logN- 
logS to those predicted by the r ecent XRB model of 
IGilli. Comastri fc Hasinger l (|2006h . This model makes 
use of the most recent observational constraints on the 
AGN populations and includes a conspicuous fraction 
of Compton thick AGN which, however, are not ex- 
pected to significantly contribute to the XMM-COSMOS 
counts. In the 0.5-2 keV band a direct comparison 
of our data with the model shows a la agreement at 
the bright end. At the faint end the model predicts 
a slightly higher normalization when compared to most 
of the plotted data, including ours. A similar behavior 
is observed in the 2-10 keV band. It is worth notic- 
ing that in the model the average unabsorbed power- 
law spectral index of the sources is < T >— 1.8 in the 
flux interval sampled by XMM-CO SMOS (see Fig. 19 in 
IGilli. Comastri fc Hasingerl (|2006D 1. Since in our data 
analysis we assumed < T > = 1.7 we expect in this band 
a slight (—10%) underestimation of the fluxes when com- 
pared to those of the model. This effect is almost negli- 
gible (i.e. <5%) in the other bands investigated here. 

By integrating our best fit 2-10 keV logN-logS be- 



tween infinite and zero, and assuming that the slope 
of the "real" logN-logS remains constant down to low 
fluxes, we estimated the total contribution of AGNs to 
the XRB. We predict a total flux of AGNs in the XRB 
of 1.25xl0" n erg cm" 2 s" 1 deg" 2 . This value is —40% 
lower than that measured bv lDe Luca fc Molendil (2004), 
and — 25% smaller than those obtained by integra- 
tion th e model logN-logS of IGilli. Comastri fc Hasingerl 
(2006) which predicts a flux of — 1.66xl0" n erg cm" 2 
s _1 deg -2 . This discrepancy between our predic ted flux 
and that of IGilli. Comastri fc Hasingerl (|2006f ). could 
arise by the fact that, in our measurement, we consider 
that all the sources have the same spectrum and from 
statistical uncertainty of the logN-logS parameters. By 
assuming an average spectral index < T > = 1.4 for all our 
sources, we obtain a value for the total flux of the AGNs 
of — 1.48x 10 -11 erg cm" 2 s" 1 deg" 2 which is consistent, 
within the statistical uncertanties, with the prediction of 
the model. Considering the total flux of the XRB pre- 
dicted by the model and our estimate from the logN-logS 
distributions, in the 2-10 keV band, XMM-COSMOS re- 
solves the —55-65% of the total flux of the XRB into 
discrete sources. 

It is interesting to observe how in the 5-10 keV band 
our data are in good agreement with the prediction of 
the model. This result is particularly important since in 
this band it is expected the major contribution of highly 
absorbed AGN, which are an important ingredient of the 
XRB models. A detailed analysis of the spectral proper- 
ties of the br ightest X-ray sources in XMM-COSMOS is 
presented bv lMainieri et al.1 (|2007l ). 

6. SAMPLE VARIANCE 

The amplitude of source counts distributio ns varies sig- 
nificantly among different su rveys (see e.g. lYang etall 
l2003t ICappelluti et al.1 l2005l and references therein). 
This "sample variance", can be explained as a combi- 
nation of Poissonian variations and effects due to th e 
clustering of sources (jPeebles I [19801 : I Yang et all l200l . 
The variance of counts in cells for sources which are an- 
gularly correlated can be obtained with: 

< (N - AfQ) 2 >= Afn + Af 2 J dOi dn 2 w(9 h2 ) (9) 

where N is the mean density of objects in the sky, £1 is 
the cell size, and ui(#i,2) is the angular two-point corre- 
lation function. The first term of Eq. [9] is the Poissonian 
variance and the second term is introduced by the large 
scale structure. In order to determine whether the differ- 
ences observed in the source counts of different surveys 
could arise from the clustering of X-ray sources, we esti- 
mated the amplitude of the fluctuations from our data, 
by producing subsamples of our survey with areas com- 
parable to those of. e.g., Chandra surveys. 
The XMM-COSMOS field and the Monte Carlo sample 
fields of Section 4 were divided in 4,9,16 and 25 square 
boxes. Making use of the 0.5-2 keV energy band data, 
we computed for each subfield, the ratio of the number 
of real sources to the number of random sources. In 
order to prevent incompleteness artifacts, we conserva- 
tively cut the limiting flux of the random and data sam- 
ple to 5xl0 -15 erg cm" 2 s" 1 . At this flux our survey is 
complete over the entire area. In order to avoid artifacts 
introduced by the missing pointings in the external part 
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TABLE 3 

Summary of the 0.5-2 keV sample variance in the 
COSMOS field. Prediction and observation at a flux 

LIMIT S;im=5xl0 -15 ERG CM -2 S _1 



FlG. 10. — The counts in cell fluctuations within the XMM- 
COSMOS field. The data are normalized to a random distributed 
field in boxes of 40', 26', 20' and 16' of side, respectively. The dashed 
lines represent the la expected fluctuation. 



of the field of view, we concentrated our analysis to the 
central 80' x 80'. In FigfTUlwe plot the ratio of the data to 
the random sample as a function of the size of the cells 
under investigation. The measured fractional standard 
deviations of the sample is reported in Table [3] Using 
Eq. [9] we computed the expected amplitu de of source 
counts fluctuations with 10(01,2) taken from lMivaii et aTl 
(2007). They computed the X -ray two-point correlation 
function in the XMM-COSMOS field and detected clus- 
tering signal with angular correlation length #0 ~1.9" 
~0.8" and -6" in the 0.5-2 keV, 2-4.5 keV and 4.5-10 
keV band, respectively. The observed slope is 7=1.8 in 
all the energy bands. 

The predicted fractional standard deviations are there- 
fore 0.13, 0.19, 0.23 and 0.28 on scales of 0.44 deg 2 , 0.19 
deg 2 , 0.11 deg 2 and 0.07 deg 2 , respectively. These val- 
ues are in good agreement with those observed in the 
sub-samples of our dataset as shown by the value of the 
fitted x 2 t° the counts in cell fluctuations (see Table [3]). 
As shown in Table [3l at this limiting flux and on the ar- 
eas considered here the main contribution to the source 
counts fluctuations is from the Poissonian noise. At the 
flux limit assumed here, the ratio a c i/a p increases from 
~0.5 on the smallest scale (16 x 16 arcmin) to ~0.85 on 
the largest scale (40 x 40 arcmin). This ratio scales as: 



eW^ccA^^-^V 3 ^)/ 2 , 



(10) 



where Af deg -2 is the surface density of the sources, Qq 
(deg) is the angular correlation length and a (deg) is 
the size of the cell. In order to estimate at which flux 
limit fluctuations introduced by the large scale structure 
are predominant, we estimate that a c i/a p would be ~1 
on the smallest scale, corresponding to a Chandra ACIS 
field of view (16 x 16 arcmin) at a surface density of 
the order of ^900 deg -2 , corresponding to a 0.5-2 keV 
flux S~ 8xl0 -16 erg cm -2 s _1 . At even fainter fluxes 
the dominant contribution to the total expected source 



Area a 


&obs h 


<T p c 




u exp 


x 2 /d.o.f. f 


arcmin 2 












40'x40' 


0.09±0.04 


0.10 


0.09 


0.13 


4.21/3 


26'x26' 


0.20±0.05 


0.15 


0.10 


0.19 


8.93/8 


20'x20' 


0.21±0.04 


0.20 


0.11 


0.23 


16.63/15 


16 'x!6' 


0.24±0.02 


0.25 


0.12 


0.28 


25.15/24 



a Size of the independent cells. 

b The observed standard deviation. 

c The predicted Poissonian standard deviation a p . 

d The predicted standard deviation due to clustering cr c i. 

e The total predicted standard deviations. 

f Value of the fitted x 2 /d.o.f. 

counts fluctuations on this area (a eX p ) comes from the 
large scale structure, therefore the contribution of sta- 
tistical fluctuations becomes less important. With the 
same procedure, we can estimate the total expected fluc- 
tuations (a exp ) and the relative importance of a p and a c i 
also for the hard band (5-10 keV), even if in this band 
we do not have enough statistics to divide our field in 
sub-samples. Using the formal best fit for 0q in this band 
found by Miyaji et al. (2007) (0 O =6" ), we find that at the 
faintest 5-10 keV flux (S~10 -14 erg cm~ 2 s _1 ) sampled 
by the XMM-COSMOS survey (M ~ 110 deg- 2 ) the ra- 
tio (Jcij u v is smaller than one on all the four scales here 
analyzed, with a total expected standard deviation of the 
fluctuations ranging from ~0.20 on the largest scale to 
^0.40 on the smallest scale. These values for a exp are 
significantly larger than those shown in Table 3 for the 
soft band, because in the hard band the surface density 
of sources is lower and the angular correlation length is 
higher than in the soft band. 

This analysis is at least qualitatively consistent with 
Figures 8 and 10, which show a significantly larger dis- 
persion in the data from different surveys in the hard 
band than in the soft band. Moreover, the results here 
discussed are also consistent with the observed fluctua- 
tions in the deep Chandra fields (see, for example, Bauer 
et al. 2004). Large area, moderately deep surveys like 
XMM-COSMOS are needed to overcome the problem of 
low counting statistics, typical of deep pencil beam sur- 
veys, and, at the same time, to provide a robust estimate 
of the effect of large scale structure on observed source 
counts. 

As a final consideration, we tried to compute the ex- 
pected intrinsic variance of XMM-COSMOS. This esti- 
mate must be made with care since we have only one 
sample on this sc ale. Assuming that the angular corre- 
lation function of iMivaii et al.l ( 20071 ) was universal, the 
residual uncertainties on the source counts are estimated 
to be <5-6% in the 0.5-2 keV and 2-10 keV bands, and 
of the order of the 10% in the 5-10 keV band. 



7. SUMMARY 

The data analysis of the first run of observations of 
the XMM-Newton COSMOS wide field survey has been 
presented. A total of 1390 point-like sources are de- 
tected on a contiguous area of about 2 deg 2 down to 
fluxes of 7.2xl0 -16 erg cm -2 s _1 , 4.0xl0 -15 erg cm -2 
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s _1 and 9.7xKT 15 erg cm" 2 s _1 in the 0.5-2 keV, 2-10 
keV and 5-10 keV energy bands, respectively. The de- 
tection procedure was tested through Monte Carlo simu- 
lations which confirmed the high level of accuracy in the 
determination of the source properties (aperture photom- 
etry and positioning) and allowed us to keep statistical 
biases under control. 

A robust estimate of X-ray source counts at both soft 
and hard energies, obtained thanks to the large num- 
ber of sources detected in the XMM-COSMOS survey, 
is presented in this paper. The differential logN-logS was 
fitted with a broken power-law model in the 0.5-2 keV 
and 2-10 keV energy bands, and with a single power-law 
model in the 5-10 keV energy band. In the soft 0.5-2 keV 
band, already extensively covered by ROSAT, XMM- 
Newton and Chandra surveys over a range of fluxes en- 
compassing those sampled by the COSMOS survey, our 
results are in excellent agreement with previous analy- 
sis (see Fig [7]), providing an independent evidence of the 
validity of our data analysis procedure. 

The large number of X-ray sources of the COSMOS 
survey allowed us to constrain with unprecedented accu- 
racy the logN-logS parameters in the 2-10 and 5-10 keV 
energy ranges over a range of fluxes which were previ- 
ously poorly constrained. Most importantly, in the hard 
5-10 keV band, we were able to fill the gap between 
the deep Chandra surveys in the CDFS and CDFN and 
shallower large area surveys. The deviations from other 
surveys, which are, however, less than 30%, have been 
explained in terms of low counting statistics of pencil 
beam surveys, and partially by the effect of large scale 
structure. The major step forward in the determina- 
tion of hard X-ray source counts achieved thanks to the 
XMM-COSMOS survey will provide an important ref- 
erence point for the study of the AGN demography and 
evolution especially with applications to obscured AGN. 
More specifically, the evolutionary properties of the ob- 
scured AGN can be tightly constrained, since they are in- 
deed very sensitive, according to the most recent model 
of the X-ray background (Gilli, Comastri & Hasinger 
2006), to the shape of the hard X-ray source counts 
around the break flux, which is precisely where the COS- 
MOS data play a key role. In this context, we compared 
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our results to the most recen t pred ictions of the model by 
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The second pass of the XMM-Newton observations in 
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expected to significantly increase the total number of X- 
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it will make possible the study of AGN evolution and 
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they reside with unprecedented detail. The COSMOS 
field has been granted 1.8 Msec observation with Chan- 
dra in its central square degree (C-COSMOS, P.I. Martin 
Elvis). The joint Chandra and XMM-Newton observa- 
tion will provide an unprecedented lab for AGN physics. 
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